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Deep multi-frequency radio imaging in the Lockman Hole: II. The 
spectral index of submillimetre galaxies 
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ABSTRACT 

We have employed the Giant Metre-wave Radio Telescope (GMRT) and the Very Large Ar- 
ray (VLA) to map the Lockman Hole. At 610 and 1,400 MHz, we reach noise levels of 15 
and 6 /iJy beam -1 , respectively, with well-matched resolutions (^5 arcsec). At this depth we 
obtained reliable detections for about half of the known submm galaxies (SMGs) in the field. 
For radio-identified SMGs, which are typically at z ~ 2, we measure a mean radio spectral 
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index of a, 



1,400 
610 



-0.75 ± 0.06 (where S v oc v a ) and standard deviation of 0.29, between 



; 1 INTRODUCTION 



SMGs were discovered in the late 1990s (Smail etal.|[l997h us- 
ing the Su bmillimetre (submm ) Common-User Bolometer Array 
(SCUBA - lHollandetai]|l999h on the 15-m James Clerk Maxwell 
Telescope (JCMT). SCUBA had been designed to exploit the so- 
called "negative A"-correction" in the submm waveband, allowing 
the detection of very distant dusty galaxies almost unbiased in red- 
shift up to 2 ~ 10. The discovery of these massive, rapidly star- 
forming galaxies (star-formation rate - SFR ~ 1O 3 M0 yr -1 ), each 
potentially capable of creating a massive elliptical galaxy within 
1 Gyr, has provided a powerful motivat ion for improving m odels of 
galaxy formation and evolution (e.g. lSwinbank et al][20 08). SMGs 
are believed by many to be the parent population of present-day 
elliptical galaxies and contribute (Sssu/jm > 3mJy) for approxi- 
mately 20 per cen t of the cosmic IR background (CIRB) at 850 /im 
(Eales et al. 1999). Understanding the nature of these galaxies has 
been extremely challenging, a situation influenced by the poor res- 
olution of submm images and the resulting confusion at noise levels 



approximate rest-frame frequencies of 1.8 and 4.2 GHz. The slope of their continuum emis- 
sion is indistinguishable from that of local star-forming galaxies and suggests that extended 
optically-thin synchrotron emission dominates the radio output of SMGs. Cooling effects by 
synchrotron emission and Inverse Compton (IC) scattering off the cosmic microwave back- 
ground (CMB) do not seem to affect their radio spectral energy distributions (SEDs). For 
those SMGs judged by Spitzer mid-infrared (-IR) colours and spectroscopy to host obscured 
active galactic nuclei (AGN), we find a clear deviation from the rest of the sample - they 
typically have steeper radio spectral indices c^q 00 ^ —1.0. These findings suggest these 
mid-IR-/AGN-selected SMGs may have an intrinsically different injection mechanism for 
relativistic particles, or they might reside in denser environments. This work provides a reli- 
able spectral template for the estimation of far-IR/radio photometric redshifts, and will enable 
accurate statistical if -corrections for the large samples of SMGs expected with SCUBA-2 and 
Herschel. 

Key words: galaxies: high-redshift — Galaxies, galaxies: active, starburst — Galaxies, radio 
continuum — Sources as a function of wavelength, submillimetre 



required to detect significant numbers of SMGs, and also by their 
i ntrinsically high redsh ifts (mean z ~ 2.2 for radio-detected SMGs 
- IChapman et all2 005) and dusty nature which make optical detec- 
tions difficult - radio and IR detections have b een usually employed 
to address these issues (e.g. llvison et alj|2007l) . 



Radio imaging has p layed a key role in characterising SMGs 
(e.g. Ilvison et aU |2002|) providing a high-resolution proxy for 
the rest-frame far-IR emission via the far-IR/radio co rrelation 
faelou et alJl 19851 : lAppleton et alj|2004l ; llbar et alj|2008l) . The ra- 
dio waveband is a relatively unexplored part of the SED where cur- 
rent studies have adopted a canonical power-law form, S v oc v a , 
based on local star-forming galaxies with a — —0.7 or —0.8 
dCondon|[T992l). and a rela t ively small dispersion, Act ~ 0.25. 
However, Hunt & Maiolind J2005h questioned this assumption for 
SMGs, arguing that the adoption of a = —0.7 is not based on 
any relevant observational evidence. Flatter radio SEDs - typi- 
cal of blue compact dwarf galaxies or some powerful ultralumi- 
nous IR galaxies (ULIRGs, e.g. Arp220, with afitcll = -0.41; 
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ISmith et"al] 1 1998h - were preferred in iHunt & Maiolincl's photo- 
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metric redshift analyses of SMGs. Indeed. IClemens et al.l (2008) 
recently measured a mean spectral index of a = —0.5 between 
1.4 and 4.8 GHz for a sample of local ULIRGs (the local 'cousins' 
of SMGs), although this steepens to a = —0.7 and —0.8 between 
4.8-8.4 and 8.4—22.5 GHz, respectively. The presence of dominant, 
compact, flat-spectrum, synchrotron-self-absorbed AGN cores in 
the SMG population cann ot be ruled out either. 

IKovacs et al. I J200d> employed 350-/im observations of 15 
radio-detected SMGs with known redshifts to constrain their char- 
acteristic temperatures. In doing so, they claimed to see the first 
signs of a deviation from the far-IR/radio correlation exhibited lo- 
cally, q L = log{Z/ F iR./([4.52 THz]Li, 400MHz )} ~ 2.3 (where 
Lfih. is the total dust luminosity and I/i,4oomhz is the rest-frame 
radio lumin osity), the ubiq uity of which has surprised astronomers 
for decades. IKovacs et alj suggested that SMGs are over-luminous 
at radio frequencies with respe ct to the far-IR , findin g a mean 
q L » 2.14 ± 0.07. Modelling bv lSwinbank et alj j20o3) suggested 
that this value is not a result of sam ple selection b ut of real evolu- 
tion of the far-IR/radio correlation. IKovacs et alj noted, however, 
that correcting the radio spectral indices by ~0.35 would move 
SMGs back onto the far-IR/radio correlation (see Jj4). 

Theoretically, at the redshifts typical of SMGs we expect 
to observe steeper radio spectral indices due to t he shift of the 
frequ ency break produced by ageing effects (e.g. ICarilli & Yunl 
Il999h. and/or IC scatte ring off the CMB - going as (1 + z) A 
(e.g. lKlameretalJl2006h . Indeed, bright radio samples selected on 
the basis of their steep (a < —1.0) spectral indices - so-called 
ultra-steep-spectrum (USS) sources - have bee n found to con- 
tain a large number of high-redshift galaxies (e.g. lDe Breuck et"al1 
2000). Given their high redshifts, their large masses and their prob- 
able relation to galaxy formation in proto-clusters, we might ex- 
pect to find similarities in radio spectral indices between SMGs 
and USS sources. We note, however, that there are major dif- 
ferences between these populations: radio-identified SMGs are 
much fainter than typical USS radio galaxies and star forma- 
tion (r ather than AGN) clearly provides the majority of their 
power dFraver et alj 1 1998]) . The fraction of SMGs with obvious 
AGN is relatively small, although Compton-thick AGN may con- 
tribute significantly to the bolometric emission in some cases. Dif- 
ferent approaches have been used to measure the actual AGN 
contribution; spectroscopic/photo metric mid-IR diagnostics (e.g. 



Menendez- Delmestre et alj l2009l) . high resolution radio imag- 



ing (e.g. iBiggs & Ivisorj l2008h . X-ray detections/stacking (e.g. 



I Alexander et alj 2003 ). and d eviations from the FIR/radio corre- 
lation (e.g. lKovacs etal]|2006l) . 

We thus have strong motives for exploring the radio spectral 
index of SMGs: to provide an independent description of their ra- 
dio SEDs, to enable more accurate radio if -corrections, to better 
test whether the far-IR/radio correlation evolves, and to find pos- 
sible new AGN diagnostics. Here, we report deep, dual-frequency, 
matched-resolution radio observations in the Lockman Hole that 
were designed to address these issues. 



2 MULTI- WAVELENGTH OBSERVATIONS IN THE 
LOCKMAN HOLE 

2.1 Parent catalogue of submm galaxies 

As one of the SCUBA Half Degree Extragalactic Survey 
(SHADES) fields, 485 arcmin 2 of the Lockman Hole was observed 



with SCUBA on JCMT at 850 fim, producing images with a FWHM 
beamsiz e of 14.8 arcsec. In this work, we use the catalogue pro- 
vided bv lCoppin etai] ( l2006h . comprising 57 sources detected with 
a signal-to-noise ratio, SNR 3.5 in more than one independent 
reduction of the data. The sample includes only sources with a low 
probability (^5 per cent) that the true, deboosted flux density is 
lower than zero (i.e. that the source is entirely spurious). 

The Lockman Hole was also observe d using the A stronom- 



ical Thermal Emission Camera (AzTEC - Wilson et alj|2008l) on 



JCMT at 1,100/xm dAustermann et al. 2009). This survey covers 
a large area: 1,115 arcmin 2 , to an r.m.s. of ^0. 9-1.3 mjy beam -1 
(~18 arcsec FWHM). In this work, we use the lAustermann et alj 
sample of 51 sources with SNR ^ 3.5 an d P(5i.ioo um ^ 
mjy) < 0.05, i.e. the same definition used by Cop pin et all 

The MAx-Planck M illimeter BOlometer array (MAMBO- 
117 iKrevsa et alj 1 19991) on the Instituto de Radioastronomia 
Milimetrica (IRA M) 30-m telescope h as also been used to image 
the Lockman Hole. iGreve et"ai". I d2004h acquired data at 1,200 /im, 
covering 197 arcmin 2 down to an r.m.s . of 0.6 mjy beam -1 
(11 arcsec, FWHM). We use the IGreve et alj catalogue which con- 
tains 23 sources with SNR ^ 3.5. 

The combination of AzTEC, SCUBA and MAMBO imaging 
makes the Lockman Hole the largest area of sky mapped at 850, 
1,100 and 1,200 /im. Combining the three samples results in an 
overall sample of 1 1 1 SMGs (20 of which were detected in more 
than one sample) which we have adopted for cross-matching with 
ou r radio catalogues. W e opted not to use the Bolocam catalogue 
bv lLaurent etliiT ( l2005l) because of the confusion resulting from the 
large beam (~31 arcsec FWHM). 



2.1.1 AGN in the sample 

Spectroscopic ana lyses using Spitzer's IR Spectrograph (IRS - 
iHouck et al]|2004l) to identify the mid-IR power sourc e in SMGs 
(e.g. Pope et al. 12008; Menendez -Delmestre et alj|2009h have con- 
firmed earlier indications that they are usually powered by ex- 
treme star formation rather than by AGN activity. Indeed, based 
on co mposite median template of a sample of 24 radio-identified 
SMGs, Menend ez-Delmestre et alj find a typical AGN contribution 
of < 32 per cent to the bolometric mid-IR emission. 

However, some of the AzTEC and SCUBA sources in the 
SHADES fields were selected for IRS observations using colour 
criteria desig ned to preferentially select highly obscured AGN 
(Copp hTet all in preparation). A sample of six such sources was 
targeted in the Lockman Hole, from which four have a dominant 
(>50-per-cent) AGN contribution to their mid-IR spectra. We use 
this sub-sample to find possible deviations in the behaviour of radio 
spectral indices (see ^3). 



2.2 GMRT and VLA imaging 

We have obtained deep radio observations in the Lockman Hole 
using the GMRlQ and the National Radio Astronomy Observa- 
tory's (NRAO) VLaB operated at 610 and 1,400 MHz, respec- 
tively. The two image mosaics have r.m.s. noise levels in the central 
regions of 15 and 6/iJy beam -1 and synthesised beams (FWHM) of 



1 GMRT is run by the National Centre for Radio Astrophysics of the Tata 
Institute of Fundamental Research. 

2 NRAO is operated by Associated Universities Inc., under a cooperative 
agreement with the National Science Fundation. 



7. 1 x 6.5 and 4.3 x 4.2 arcsec 2 , respectively (see a full description 
in llbar et al]|2009l) . 

In this work, we investigate the radio spectral indices of the 
SMG sample described in £|2.1| exploiting our Lockman Hole ra- 
dio data down to a peak-to-noise ratio, PNR ^ 3 - where peak 
and noise refer to the maximum value of a 2-D Gaussian fit and the 
local r.m.s., respectively. Given the resolution of our radio imag- 
ing (4-7 arcsec, FWH M) and the typical kpc scale of SMGs (e.g. 
Biggs & Ivisonl i2008). we have considered the SMG counterparts 
to be unresolved. We have extracted flux densities with the source 
size in the SAD Gaussian fitting routine within the Astronomical 
Image Processing Syste rr0 (MVS) set to the size of the synthe- 
sised beam (DOWIDTH = —1), which results in a less uncertain in- 
tegrated flux density estimate because of the lower number of free 
parameters in the Gaussian fit. Monte Carlo simulations based on 
fake Gaussians injected in a residual map (with no > 5-cr peaks) 
using IMMOD and extracted using SAD show a clear improvement 
(~1.2x less scatter in the out:in flux density ratios at PNR < 5) 
by fixing the synthesised beam in this way. Fewer outliers are seen 
as well. 

There is an important caveat concerning the use of these faint, 
low-PNR sources: a very large fraction of the sources detected at 
faint flux densities are due to the background noise. To determine 
the prevalence of spurious sources in the samples, we inspected 
the number of 3-<x sources extracted from the inverted residual 
map. The cumulative number of negative sources, A r ncg (>S), in- 
creased rapidly toward faint flux densities, overwhelming the num- 



ber counts of the 'positive' sources at Si, 4 



40 /iJy and 



S6ioMHz & 120 /iJy. This indicates that an extrapolation of the 
number counts to Si, 400MHz £ 100/iJy (S6iomhz ;£ 300 /uJy), 
based on sources with higher flux densities, will thus underesti- 
mate the number of extracted sources in the field. This is relevant 
as we need to know the probability of random associations between 
SMGs and faint radio sources. 



2.3 Radio identifications 

Associating SMGs with radio counterparts has usually been a ccom- 
plished using probabilistic arguments (e.g. llvison et alj|2007l) . We 
have chosen to cross-match the parent SMG catalogue with our ra- 
dio maps ( i]2.2t using a P < 0.05 selection criterion, meaning 
that the probability o f a spurious association is less than 5 per cent 
dPownes etaiTl 986). This probability has two dependencies - the 
radius within which we search for radio counterparts and the cu- 
mulative radio source number count distribution per surface unit, 
N(>S). The search radius was chosen to be 3 x the uncertainty in 
the submm position (in Right Ascension o r Declination), as esti- 
mated using Eq. B22 of llvison et 10 J2007t) , i.e. dependent on the 
submm SNR, with 12.0, 14.4 and 8.8 arcsec for SNR = 4 SCUBA, 
AzTEC and MAMBO detections; we expect 1 per cent of the true 
counterparts to lie outside our search area. For the number count 
dependency, we used all ^ 3-cr radio detections to find the true 
probability of random association. Nevertheless, if the source was 
detected at both radio wavelengths (unlikely to be noise) we used a 
negative-subtracted cumulative radio source number count distribu- 
tion instead, N Buh (>S) = N(>S) - N ncg (>S), in order to give a 
better proxy for the number of real sources. We note, however, that 
the power of this P-identification procedure is lessened for very 
large submm beams and/or low SNRs because the purely statistical 



see http://www.aips.nrao.edu/cook.html for a description of tasks. 



approach is unable to discriminate against unrelated radio-bright 
sources (preferentially radio-loud AGN) at large separations. 

We adopted a conservative approach, eliminating from the 
analysis all SMGs with more than two radio sources that could af- 
fect the measurement of the radio spectral index, thus biasing our 
sample against close mergers. A final inspection by eye was used 
to eliminate sources with confused associations. The final sample 
comprised 44 SMGs (from a total of 111) with single and reliable 
radio counterparts in at least one radio waveband (41 and 37 detec- 
tions by the VLA and GMRT, respectively). 



3 THE RADIO SPECTRAL INDEX OF SMGS 
3.1 Observations 

In Fig. [T] we show the radio spectral indices of the 44 radio- 
selected SMGs as a function of flux density at 1,400 MHz. The 
vast majority of these sources were detected at low flux densi- 
ties, Si, 4oomhz ^ 150 /iJy, close to the detection limits of even 
the deepest radio imaging. 

The study of the radio spectral index has a clear bias intro- 
duced by the detection threshold in each map. The solid black line 
in Fig. Q] shows the restriction for detecting flat-spectrum sources 
at faint 1,400-MHz radio fluxes, based on a minimum 3-<r source 
at S6iomhz ~ 45 /j,Jy in the deepest region of the 610 MHz image 
(note that noise varies within a factor of < 2 with position in the 
maps). In order to reduce the effect of this bias, we have restricted 
the statistical analysis to those sources with Si,4oomhz > 40 /J,]y 
only (see shaded area in Fig.[TJ. From this sub-sample, 28 were de- 
tected at both radio wavelengths, while one and three were detected 
only by the VLA and GMRT, respectively. We find - for those 
SMGs detected at both radio frequencies - a bootstrapped median 
radio spectral index of —0.71 ± 0.04, a mean of —0.75 ± 0.06 
and standard deviation of 0.29. In the data, there is no clear corre- 
lation between o^g 00 , radio and submm flux densities, suggesting 
the radio spectral indices of the radio-identified SMGs shown in 
this work may provide a good template for SMGs below the cur- 
rent submm confusion limit. 

By looking at the sample of four AGN s selected via mid-IR 
spectroscopic analysis ( ij2.1.11 Coppi rTet all in preparation), we 
find all of them have a reliable a^g 00 measurement ( i)2.3t . Al- 
though a small sample, restricted to bright mid-IR fluxes, their 
radio properties clearly deviate from those of the other SMGs 
(filled black circles in Fig. [TJ. One is a flat-spectrum source, 
Q 6K>° = —0.05, while the other three present steep radio spec- 

1,400 , -1 n 

tra, a 6 ' 10 ^ -1.0. 

The different nature of these sources is also evidenced by the 
detectio n of the steepest radi o source (mid-IR AGN-dominated) in 
X-rays dBrunner et alj|2008h showin g clear eviden ce of a highly 
obscured AGN - in agreement with ICoppin et Us criteria - we 
note this is the only SMG with a clear X-ray detection. Flat radio 
spectra are also related to blazar or powerful compact AGNs and 
are usually interpreted as the synchrotron self-absorbed base of a 
"young" expanding jet that fuels larger scale radio emission lobes. 



4 DISCUSSION 



The observed radio spectral indices for SMGs, at approximately 
rest-frame 1.8 - 4.2 GHz for z ~ 2, are similar to either those 
found in local star-forming galaxies ( Condon 1992 ) and to the bulk 
of the sub-mjy radio sources found bv llbar et al.U2009h . The flux 




Figure 1. The radio spectral index, between 610 and 1,400 MHz, of a sample of 44 SMGs selected from the Lockman Hole as a function of observed radio 
flux density at 1,400 MHz. Plus, cross and square symbols represent unique AzTEC, SCUBA and MAMBO sources (see § 12.3) . Error bars are estimated from 
the uncertainty in flux density at each frequency. Up per and lower limits are 3 a thresholds. Large filled circles are SMGs dominated by AGN activity as 
judged by Spitzer mid-IR colours and spectroscopy I Cop pinet all in preparation). A cross shows the source is detected in X-rays. The solid line denotes 
the threshold in radio spectral index based on a 3-ct source with SeiOMHz = 45 /jjy. The mean and standard deviation (using only Si^ooMHz > 40 fiJy 
sources) are shown in dark and light grey, respectively. The box at the bottom-right corner shows the histogram of radio spectral indexes for those sources with 
Si,400MHz > 40 /Jy; solid, dashed and dotted lines are detections, upper limits and the normal distribution fit, respectively. 



density range at which the vast majority of the SMGs are found, 
Si,4oomhz < 150 /iJy, suggests they may compose a significant 
fraction of the radio population responsible for the flattening of the 
Euclidean radio number counts. 

Assuming the physics of the synchrotron emission does not 
evolve as a function of redshift, these results suggest that non- 
thermal optically thin synchrotron emission dominates in the ma- 
jority of SMGs. This sug gests magnetic fields a re large enough 
to retain cosmic rays (see Thompson et al. 2006) and contradicts 
the idea that flat ter radio SEDs - such as those thermally dom- 
inated galaxies jHunt & Maiolind I2005T) or in synchrotron self- 
absorbed sources - are more appropriate SEDs for these massive 
star-forming galaxies. Indeed, the lack of flat spectrum sources 
suggests that their emission is mostly dominated by extended 
rather than compact radio structures, in a greement with previ- 
ous SMG high-resolut ion radio observations iChapman et al . 2004; 
Bigas & Ivison 2008) and the low mid-IR extinctions found by 
Menendez-Delmestre et al. (2009). 

The lack of evidence for spectral steepening with respect to 
local star-forming galaxies suggests that i^-correction effects in 
curved SEDs and IC scatte ring enhanced by a denser CMB radi- 
ation field at high redshift jKlamer et al.l2006r . lMilev & De Breuckl 
2008) do not have a strong influence on the observed radio spectral 
index of these high-redshift galaxies. The possibility of intrinsic 
redshift evolution to flatter radio spectra, where curved SEDs (due 
to synchrotron and IC cooling) steepen them back to the observed 
values, cannot be ruled out. 

The steeper radio spectral indices found for the mid-IR AGN- 
dominated sources are intriguing. Recent studies of USS samples 
(also at high redshift and mostly AGN) using multi-frequency radio 
observations have found no spectral curvature at higher frequencies 



in these samples. The power-law behaviour suggests steep radio 
spectr a in high-redshift A GN are intrinsic, or a product of environ- 
ment ( Bryant et alj|2009l) . The former relates to a different mech- 
anism for the injection of the initial electron e nergy distribution 
(e.g. IJaffe & Perolalll973l ; ICarilli & Yun| [l999) which determine 
the reservoir of high-energy particles - we expect to observe steeper 
radio spectra if there is no continuous injection or re-acceleration 
of relativistic particles across the galaxy. On the other hand, if the 
source is embedded in a dense medium, some of the energy of the 
particles would be released into the environ ment and stronger mag - 
netic fields may increase the cooling rate dThompson et al.l 2006). 
i.e. the radio SED would become steeper. 

Based on the fact that we do not observe steeper radio spec- 
tral indices in the bulk of the SMGs, these findings may suggest 
that the synchrotron emission in AGN is fundamentally different 
to the emission seen in massive star-forming galaxies. If we as- 
sume the a vs red shift correlation found in brighter radio samples 
jDe Breuck et alj|2000h , it is possible that the synchrotron radia- 
tion emitted by nuclear activity evolves differently as a function of 
cosmic time than that emitted due to star formation. Nevertheless, 
there are major differences between USS and SMGs as pointed out 
in E 

We have checked the possibility for an overestimation of th e 
radio luminosities in the sample analysed bv lKovacs etail d2006t) . 
From the 9 SMGs in the Lockman Hole (Table 1 from their pa- 
per), 5 of them have got a reliable estimate of radio spectral index, 
Q 6K> 00 ~ —0.68 ± 0.26 (bootstrapped median). Based on these 
detections, it seems that their sample is not particularly affected by 
peculiar radio /^-corrections that may explain the observed devia- 
tion from the local far-IR/radio correlation. 

A future analysis - using radio observations from the LOw 



Frequency ARray (LOFAR) and Expanded VLA to cover a wide 
range of frequency, for samples with well-constrained redshifts - 
will allow us to give a much better description of the synchrotron 
emission processes in massive star-forming galaxies and radio faint 
AGN. An independent assessment can also be obtained using high- 
resolution radio observations in order to determine the physical ex- 
tent, morphology, luminosity and brightness temperature of their 
emitting regions. 



5 CONCLUSIONS 

We have analysed the radio spectral indices, agi^ 00 , based on well- 
matched observations with GMRT at 610 MHz and the VLA at 
1,400 MHz, of a sample of SMGs selected in the Lockman Hole 
using the AzTEC, SCUBA, and MAMBO cameras. 

We have created a sample of 44 SMGs (from a total of 1 1 1 
unique SMGs) which have secure radio identifications and reliable 
estimates of radio spectral index. We report how deep GMRT ob- 
servations have resulted in a detection rate similar to that usually 
accomplished with the VLA. We find a mean value of c^q 00 = 
—0.75 ± 0.06 (see Fig. [TJ, suggesting the majority of SMGs are 
dominated by optically thin synchrotron emission from extended 
radio structures; star-forming regions and/or AGN lobes. The dis- 
tribution of spectral indices for SMGs is indistinguishable from that 
of local star-forming galaxies. This supports the idea that SMGs 
are dominated by processes typical of dusty star-forming galaxies 
rather than those seen in AGN. If we assume that there is no in- 
trinsic redshift evolution for the synchrotron emission mechanism, 
we find that ageing effects and IC scattering off the CMB does not 
seem to largely affect - steepen - the bulk of the observed radio 
spectral indices. We do not observe any dependency of the radio 
spectral index on the submm/radio flux density ratio, or on submm 
or radio flux densities, which suggests the radio spectral properties 
found in this work should also be representative of the SMGs below 
the current confusion limits. 

We have found that the radio spectral indices of SMGs with 
an AGN-dominated mid-IR spectrum deviate from the bulk of the 
sample. They typically have either steeper radio spectral indices, 
a 6io° & -1-0 -similar to USS sources found in bright radio sam- 
ples - or flat spectra. These findings suggest a different mechanism 
responsible for the synchrotron emission in these sources: either 
related to the injection mechanism of relativistic plasma, or a prod- 
uct of the environment. Their deviation to steeper spectral indices 
could be an evidence for a different cosmic radio evolution with re- 
spect to massive star-forming galaxies. Besides of the small sample 
of sources, radio spectral index thus represents a useful selection 
criterion for identifying AGN in high-redshift SMGs. 

Looking to the future, these results can be used for statistical 
analyses of the large SMG samples anticipated from SCUBA-2 and 
from Herschel. Our parameterisation of the radio SED of SMGs 
will allow reliable if -corrections for flux densities in the ~l-GHz 
regime, and will alleviate degeneracies faced by long-wavelength 
photometric redshift techniques. 
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